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I. INTRODUCTION 



Einstein's classical theory of gravitation has been successfully applied to the solution 
of problems concerning astronomical scales of length and time, such as the formulation of 
cosmological models. In the development of cosmology much effort is being focused on 
models in which the universe was expanding at very early stages of its evolution. There 
is observational evidence that at present the Universe seems homogeneous and has been 
highly isotropic since the recombination era. In fact, observations of the cosmic microwave 
background radiation, galaxies, and other astronomical objects reveal that our universe, on 
a very large scale, is remarkably uniform and is currently under accelerated expan- 

sion It follows that on large scales the accelerated expansion is the mechanism to 

solve the problem of the apparent large-scale isotropy of the universe and its behaviour is 
approximately described by a Friedmann-Robertson- Walker (FRW) model. 

The claimed isotropy of the universe is deduced from reliable measurements, which are 
the most accurate ones in cosmology. However, these results have been established for times 
after the era at which the universe became transparent to radiation. Their extrapolation to 
earlier times and, in particular, near the Planck or string scale is totally unfounded. There 
are theoretical ar gum ents that sustain the existence of an anisotropic phase that approaches 
an isotropic one [lOj. In addition, the desire to avoid the assumption of special initial 
conditions implied by the FRW models suggests the study of more general cosmological 
models such as anisotropic and spatially homogeneous ones. 

It is generally accepted that the early universe was characterized by a highly irregular 
expansion mechanism. Therefore, it is interesting to study more appropriate cosmological 
models in which anisotropics, existing at early stage of the expansion, are damped out in 
the course of the evolution. This investigation has increased since in Ref. ^lj it was shown 
that the creation of scalar particles can dissipate the anisotropy as the Universe expands. 
A physically acceptable cosmology should provide a mechanism to achieve an accelerating 
expansion at the present time. 

A Bianchi type-I universe is of particular interest because it is one of the simplest models 
that contains special isotropic cases and allows arbitrarily small anisotropic levels [3]. Since 
the actual Universe is surprisingly isotropic, it makes a suitable candidate for studying the 
possible effects of an anisotropy in the early Universe on present-day observations. 



2 



The solutions of the low-energy string cosmological effective action are by their nature 
anisotropic and there is evidence that the dynamics of the early universe may have 
been profou ndly influenced by the presence of spatial anisotropies just below the Planck or 
string scale jl4j |. In a new family of exact spherically symmetric solutions was found in 
the model with a one-component anisotropic fluid for a general equation of state. In this 
context, it is important to analyze the role of different sources of anisotropy in the Einstein 
equations. 

Internal symmetry has proven to be a useful tool for field theory in flat spacetime to un- 
derstand charge conservation, isospin, etc., and gauge theories have played an important role 
in understanding high-energy phenomena. Any theory that attempts to unify the different 
physical interactions must discuss the two following points: the structure of the dynamical 
interactions among the fundamental fields, and the nature of the internal symmetry group 
(ISG) that governs these interactions. We are interested in the second point, not with a 
phenomenological motivation, but for a more fundamental reason, e.g., the properties of 
the Einstein gravitational action itself. As far as we know, few efforts have been made to 
investigate this kind of symmetry in the context of general relativity. 

In |l6[ we exploited and investigated an unusual kind of internal symmetry of the Einstein 
field equations in homogeneous spacetime, e.g., the form invariance symmetry, which relates 
the expansion rate (geometrical variable) with the energy density and pressure of the perfect 
fluid (source variables). It has proved to be a useful concept to give a satisfactory explanation 
of the assisted inflation mechanism, appealing to the cooperative effect of adding the energy 
density into the Friedmann equation. 

Below, we investigate the internal symmetries of the Einstein equations for the particular 
case of an expanding universe described by anisotropic Bianchi type-I metrics throughout 
its whole evolution. 

The paper is organized as follows. In Sec. II we find the linear transformation of the 
expansion rates which preserves the Einstein gravitational action for the Bianchi type-I 
spacetime and present the three-dimensional internal space. In Sec. Ill we show that a vector 
representation of the ISG exists where two scale factors of the Bianchi type-I cosmology with 
perfect fluid become equal. We find the general solution for a perfect fluid with a constant 
baryotropic index and investigate its asymptotic behavior. In Sec. IV exact solutions are 
found for three types of anisotropic fluids. Finally, the conclusions are stated in Sec. V. 
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II. INTERNAL SYMMETRY GROUP 



The spatially homogeneous anisotropic Bianchi type-I spacetime is described by the line 
element 



ds 2 = -dt 2 + aiitfdx 2 + a 2 {t) 2 dy 2 + a 3 (t) 2 dz 2 . 
We write the classical Einstein gravitational action Se for the metric (|T|). 



(1) 



s E [g] = J V^[h 1 + h 2 + h 3 + Hi + H 2 + H 2 + 

H X H 2 + H X H 3 + H 2 H 3 ] d\ (2) 

where g denotes the metric determinant, Hi = hi/ai are the expansion rates along the 
three spatial directions, and i = 1,2,3. Throughout this paper, we use units such that 
c = 8irG = 1. This gravitational action is invariant under a linear transformation of the 
expansion rates that preserves the following quantities: 



Q 1 = H 1 + H 2 + H 3 



(3) 



Q 2 = Hi + H 2 2 + Hj. (4) 

These invariants can be expressed in terms of the expansion scalar 9 and the shear scalar a 
as Q\ = 9 and Q 2 = a 2 + 9 2 /3. Introducing an internal three-dimensional cartesian space, we 
can see that the points of the circumference defined by the intersection of the plane © with 
the surface of the sphere (J3J) are equivalent. Hence, the linear transformations H[ = H'^Hj), 
which transform points of this circumference into points of the same circumference, are 
internal symmetries of the Bianchi type-I geometrical action. A matrix representation of 
the ISG, isomorphic to the rotation group on a plane, is given by 



M = - 

3 



( l + 2cos0 1 — cos (f) — sin <\> 1 — cos (j) + y/3 sin ^ 

1 — cos + \/3 sin l + 2cos0 1 — cos — sin , (5) 
i 1 — cos0 — sin0 1 — cos0 + a/3 sin l + 2cos0 j 

where is the rotation angle. The elements of the rotation matrix M have the additional 
properties J2i=i ~M ik = 1 and J2l=i ^ifc = 1 to preserve the quantity Qi defined in Eq. (J3j). 



The variational problem in the gravitational action for a stress tensor representing a 
fluid with energy density p and pressures pi, leads to the Einstein equations 



3# 2 = P+V, 



(6) 



- 2H = p + p + (J 2 , 



(7) 



p + 3H(p + p) + a ■ E = 0, 



(8) 



a • <? + 3# a 2 



a ■ E. 



(9) 



where the average expansion rate H and the shear vector a are the trace and the traceless 
part of the extrinsic curvature Kij of the homogeneous time slice. The former can be 



components o,i = H^ — H . The isotropic pressure is p = (pi +P2 +Ps)/3 and the components 
of the vector S are the transverse pressures Ej = pi — p. In addition, a and E transform as 
a vector in the internal space. Note that decomposing the three-dimensional internal space 
as a direct sum of one-dimensional space containing the rotation axis and a two-dimensional 
orthogonal space, we can see that H and p are invariants under the action of the ISG, while 
a and E are vectors contained in the two-dimensional orthogonal space. 
From the dynamical equation of the shear @ we get 



where the three integration constants <7jo transform as the components of a vector in the 
internal space and satisfy the transverse condition dio+^o+^so — with Cio+clo+^io = a o- 
Defining the dimensionless constants qi = u^q/aQ such that 



expressed in terms of the average scale factor a 



(— (7) 1 / 6 , as H = a/ a, and the latter has 




(10) 



qi + <?2 + <?3 
along with the functions m and pi, 



= 0, q\ + q\ + ql = q\ 



(11) 



where a = (cro/*?) 1 ^ 3 ) then the scale factors along the three spatial directions are obtained 
by integrating the shear vector (fTUj) 

di = a i0 a^im q \ (13) 

where the integration constants and the functions /ij satisfy the constraint 
010020030/^1/^2/^3 — 1- 

Due to the vector nature of the expansion rates Hi the three scale factors a\ give rise to 
the following vector representation of the ISG 

\na[ = M ik \na kl (14) 

and the internal symmetry transformations map solutions of the Einstein equations ©- 
(jHJ) into solutions of the same equations. In the next section we will use Eq. ()14j) to 
reduce the internal degrees of freedom of the Bianchi type-I metric. On the other hand, 
the transformation rule (|14j) extends the results of the previous paper [ifj], where the form 
invariance symmetry of the Einstein equations for a Friedmann-Robertson- Walker spacetime 
was shown . 



III. BIANCHI TYPE-I COSMOLOGY WITH A PERFECT FLUID 



A wide class of cosmological models have the initial stage similar to that of the vacuum 
Bianchi type-I model first considered by Kasner. Investigations of such spacetime filled with 
a perfect fluid were carried out in Q-Q- Later on in Q,Q it was shown that the 
evolution of the model with a stiff fluid is significantly different from that of models with 
dust or radiation, and the initial vacuum assumption can be violated. 

In this section we use the internal space to reexamine the Bianchi type-I cosmology with 
a perfect fluid characterized by an isotropic pressure p. For this fluid the transverse pressure 
S vanishes, so Eq. (|12J) gives \i{ =const and the shear (|T0|) turns into 

* = 5- (15) 

Then a 2 = Cq/o 6 is the contribution of the shear to Eq. (JfjJ) and the Einstein equations 
(J0J)-Q look as if they were the Friedmann-Robertson- Walker equations for two separately 
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conserved perfect fluids, one of which behaves as if it were a stiff fluid. In the g^-parameter 
representation the scale factors (|T3j) read 



Oj = a i0 a m q \ (16) 

where the three integration constants a i0 satisfy the condition aio(i2o a 30 = 1- 

As Eqs. ©-© contain an explicit dependence on the constant q 2 and each scale factor 
(fT7)|) is generated by the constants qi, two different solutions with labels and q[ satisfying 
Eq. pip can be related by an internal symmetry transformation. The linear transformation 
q\ = Mjfcgfc induces the transformation ()14|) between two different solutions of the Einstein 
equations with the same source and shear scalar. 

The ISG can be used to reduce the internal degrees of freedom of the Bianchi type-I 
metric. In fact, the scale factors generate a vector representation of this group; hence, 
for a given set of constants qi a system of coordinates exists in the internal space where two 
components of the transformed constants q[ become equal. This job is done by the matrix 



M (<?i=^) = 3 



' l-V&q 3 /q l-y/6q 1 /q l-Vdq 2 /q 

1 - VQq 2 /q 1 - vW<? 1 - VQqi/q • (17) 
\ 1 - \/6qi/q 1 - V6q 2 /q 1 - V^Qz/q ) 

Using this matrix in the transformation rule (pPfj) . we obtain the scale factors in the new 
representation 

a[ = a 10 am q/V ^, a' 2 = a 2Q am q/V ^ , a' 3 = a 30 am-^ q . (18) 

This shows that any Bianchi type-I cosmology with a perfect fluid can be mapped by means 
of an internal symmetry transformation to a locally rotationally symmetric model with the 
same source and shear scalar. To see this result explicitly we fix the constant q 2 = 2/3 in 
the first Eq. (JT2J); then the scale factors (fT7)|) become 

m Sl m S2 m Si 

a l = a 10 a ■ i /q , 0-2 = 0-20(10 . i in i a 3 = ^30^0 ■ i /q , (19) 

where the parameters Sj = qi + 1/3 satisfy the Kasner constraints 

si + s 2 + s 3 = 1, sj + s\ + s\ = 1. (20) 
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A one-parameter representation of the Kasner exponents Sj is given by 



1 



1 + 



2b 



52 = 3 



1 + 



3-6 



S3= 3 



3 + 6 



21) 



where b is the parameter. Hence, for the Bianchi type-I spacetime filled with a perfect fluid, 
we conclude that the set of solutions generated by the Kasner exponents Si are related 
to each other by internal symmetry transformations. In particular, for b = 1 we reobtain 
the locally rotationally symmetric solution (|18|) . 



General solution 



Now, we investigate a cosmological model with perfect fluid. The energy-momentum 
tensor of the fluid is T ik = piiiUu + ph ik) where p is the energy density, p the equilibrium 
pressure, u l the four- velocity of the fluid, and h%k the projection tensor h ik = g%k + UiU k . 
The components of the fluid four-velocity are u l = (1,0,0,0). An equation of state of 
the form p — (7 — l)p with a constant baryotropic index 7 is also assumed. Integrating 
the conservation Eq. (|SJ) we get the energy density p = A/a 37 , where A is an integration 



constant. In this case, introducing the variable v 
become 



-g in Eqs. © and JEJ), they 



v' 2 = 1 + 3\v 2 — = -, 

m v 

where the prime indicates differentiation with respect to the dimensionless time T 

3/2 <7ot and A = 2A/3o"q. The general solutions of both equations in 
follows. 



(22) 



are given as 



Case A > (positive energy density) 



3A 



sinh 2 (r — r ) 



1/(2-7) 



m 



tanh 



1 2/(2-7) 



(23) 



T 



2-7 



(3A) 



1/(7-2) 



[sinh r 



7/(2-7) 



dr . 



(24) 



where r is an integration constant. From Eqs. (|TUjl. and the scale factors are 
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3A 



1/3(2-7) 



cosh 



(r - to) 



4/3(2-7) 



tanh ■ 



7-0 



2^/(2-7) 



(25) 



If the baryotropic index is restricted to the range < 7 < 2, then from Eq. ([21$ we see that T 
and r have the same asymptotic limits. This allows us to investigate the evolution of the scale 
factors in two asymptotic regimes using the exact solutions (j23| ) -(|25j l . In the first regime, 



a 3 < (3\) 1 ^' y ~ 2 \ the perfect fluid is dynamically unimportant, the approximate average 
scale factor (|23|) becomes a 3 oc T, and m oc T. Hence, qualitatively the behavior of the scale 
factors (|25jl is similar to that of the vacuum Kasner solution. Then, in first approximation, 
the fluid can be considered as a test matter on a highly anisotropic background. In the 
second regime, which starts from some characteristic time where a 3 > (3A) 1 /*- 7-2 ), the fluid 
becomes dominant. In this case, the approximate average scale factor ()23|) and the scale 
factors (|2l)|) are a oc Oj oc T 2 / 37 . This result shows that, owing to the spatial isotropy of 
the stress-energy tensor, the anisotropic Bianchi type-I model evolves into a Friedmann- 
Robert son- Walker cosmology and the initial anisotropy of this simple cosmological model is 
dissipated as the Universe expands. 



Case A < (negative energy density) 



"3A Sm (r - Toj 



1/(2-7) 



m 



tan 



2/(2-7) 



(26) 



T 



2-7 



-3\) l ?^> I [sin tF /(2 " 7) dr. 



From Eqs. IfTHjl . ((2EJ) and (|2Tjl the scale factors are 



— a i0 



4 v 1/3(2-7) 



3A 



cos ■ 



4/3(2-7) 



tan 



[T ~ To) 



(27) 



2s7(2-7) 



(28) 



they represent a universe having a finite time span. 

For the particular cases A = and 7 = 2,1,0, the explicit general solutions are given 
below. 



1. X — (vacuum), 
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a 3 = T, m = T, a { = a i0 T s \ 



(29) 



,2. 7 = 2 (stiff fluid), 



a 3 = V1 + 3AT, 



m 



rpl/VT+3\ 



1 (si-!/3) 

a, = a, (l + 3A) 1/6 T3 + ^TT5r, 



(30) 



5.7 = 1 (dust), 



T 



m = 



T + 



3A 



/3A\ 



1/3 



T + 



3A 



4 \ 2/3 

T+ — 

3A/ 



(31) 



^•7 = (positive cosmological constant), 



1 v^3A^ 

^sinhVzXT, m = tanh T, 

^A 2 



1/3 



3A, 



, v3A_ 
tanh T 



cosh 2 / 3 ^ T. 
2 



(32) 



5. 7 = (negative cosmological constant), 



1 / v 73 ^! 

smv— 3 AT, m = tan — - — T, 



v^3A 



1/3 r 



a i — a i0 



,v /z 3A / 



V^SA^ 
tan — - — T 



2/3 V— 3A 
cos ' — - — T. 



(33) 



Explicit solutions can also be found for other 7 values, but for physical reasons we do not 
pay attention to them. 
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IV. ANISOTROPIC FLUID 



The kinetic and potential energy of the gravitational field that may be attributed to the 
presence of anisotropy could reach the same order of magnitude or even exceed the energy 
attributed to any other source present in the effective action at some time in the early 
evolution. Hence, the dynamics of the universe at that time could be profoundly distorted 
if one takes into account that contribution. To investigate these effects it is convenient to 
decompose the vector pressure in the internal space as a sum of a vector along the rotation 
axis vector n = (1, 1, 1) and a vector contained in a plane perpendicular to this axis, then 
V = pn + E. 

Now we can find exact solutions of the Einstein equations ©-© for three different 
sources, an anisotropic fluid with constant stress-energy tensor, an anisotropic fluid charac- 
terized by a vector pressure whose components are proportional to the energy density, and 
an anisotropic stiff fluid. 



A. Constant stress-energy tensor 

This fluid has a constant energy density p and a constant vector pressure V = (pi,P2,P3)- 
From Eqs. ©, © an d (fTUj) . we have 

{a 3 )" - (3a 3 = 0, P=\{p-p), (34) 

and 

°i = + (35) 
a 6 p 

where Sj = Pi — p are the components of the transverse pressure. On the other hand, the 
conservation Eq. (jSJ) gives 

" + P+ M§ + J = °< 22 = 51 < 36) 

Choosing the integration constants so that the vectors a® and S are perpendicular, the third 
term in Eq. (|36|) vanishes. This choice, which is invariant under the action of the ISG, 
selects solutions for which the equation of state is 
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or in terms of the pressures pi 



P 2 = P 2 + 



(37) 



2 (pip 2 + P1P2 + P2P3) =P 2 i+pI+pI + 3p 2 . 



(38) 



When the energy density p is greater than the transverse pressure £, we expand the equation 
of state ()37|) in powers of £/p and p~ ±(p+£ 2 /3p +...). The negative branch could describe 
a modified Chaplygin gas (see |20|); by a Chaplygin gas we mean a perfect fluid with the 
equation of state p = —A/p where A is a constant and p < 0. This equation of state has 



raised a certain interest recently due to its many interesting features 



2l|. It has a connection 



with string theory and can be obtained from the Nambu-Goto action for (i-branes moving 



in a (d + 2)-dimensional spacetime in the light-cone parametrization 22| . In addition, the 
Cha plyg in gas is the only fluid which, up to now, admits a supersymmetric generalization 
24j . In the above approximation the constant A = S 2 /3 would represent the first 



contribution of the anisotropy pressure to the isotropic average pressure. 
From Eqs. (|T2J) and (|3*Hj) we obtain the scale factors 



a i0 m qx a $ , 



(39) 



where the variable v = a 3 satisfies the equation 



3 2 P~P 

4 p 



1 + -V 



(40) 



which has been obtained by inserting Eq. (J35)) in Eq. For p = — p = A Eq. (|4U|) 

reduces to Eq. (j22j) for 7 = 0, Eq. (|3*7|) leads to a vanishing transverse pressure, and the 
fluid becomes a perfect fluid. In other cases, introducing the parameter 3A = 2p/a^ and the 



dimensionless time T = ao^/3(p — p)/p t/2 in Eqs. (JT2J) and (|4~0|). they turn into Eq. 
for 7 = after fixing q - 



'2p//3. Inserting its solutions in the scale factors (j%2)) we get 
three different types of solutions 



1. p > and p < — p 
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di — aw 



tanh 1 
2 



: sinh \T(3 1 ] 



2p 



(41) 



2. p < and p > —p 



a 3 = — , °_ sin J— (3 1, m 



tan — — t 



di (Xi 



tan - ' — t 



sin d —p t 



(42) 



3. p < and p < p 



cosh vr/3t, m = exp 



2 tan e 



a* = a i0 exp 



2g i tan~ 1 eV / ^* 



O"0 



cosh v/ p t 



1 4.^1 
3^ 



(43) 



Solutions (|41|) and (|42|) could be considered as a generalization of the solutions (J32j) and 
(I33|) with constant cosmo logical term. In fact, when the transverse pressure vanishes both 
sets of solutions are the same. 



B. pi oc p 

Here we will find exact solutions for a cosmological Bianchi type-I model with an 
anisotropic fluid characterized by a vector pressure whose components are proportional to 
the energy density and investigate their asymptotic behavior. The equation of state is 

V={Q- n)p, (44) 
with Q = (71,72,73) a constant vector index that can be written as 

Q = Q n + Q* = 1 n + f, (45) 
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where < 7 < 2 is the average index and f is the transverse index. The particular case of 
7 = 2, anisotropic stiff fluid, will be investigated in the next subsection. Comparing Eqs. 
and with V = pn + E, the average pressure and the transverse pressure are 



p = (7 - l)p, S = I> (46) 

For vanishing transverse index the average pressure p reduces to the isotropic pressure and 
the equation of state becomes that of a perfect fluid. 

Combining Eqs. ©-© with the first one in Eq. (j4^|) . we obtain 



Using Eqs. (fT0|) . fT2|) . (}46|) and (J4?j) . we calculate the shear vector 



ao 2T 

a 3 2 - 7 

and the corresponding scale factors after integrating Eq 



cii = a i0 m q ' a + 1-~> . (49) 

With the aid of Eq. (fT2"|) . the transverse pressure (J4T)|) . and the shear vector ()48|) we can 
integrate the conservation Eq. (JSJ) to obtain the energy density of the anisotropic fluid 

Am CT o 

P = -^jET> (50) 

where A is an integration constant and T 2 = 7i + 72 + 73 — 37 2 . Defining the effective index 
7 e = 7 + 3(2-7) an d choosing the vector cr perpendicular to f , the energy density of the fluid 
takes the form p = A/a 3> and the Einstein equation (jBJ) for the variable v = a 3 becomes 



• 2 2 



3 2 2 - 7 



1 + 

^0 



2--2- 7 . (51) 
Finally, the scale factors are obtained after inserting the solutions of Eqs. ()12j) and (JoTj) 
in Eq. (|49j). For fluids with average index < 7 < 2, the contribution of the transverse 
pressure increases the effective index. Consequently, the anisotropic pressure decelerates the 
average expansion of the universe. 
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The use of the dimensionless time T = a J3(2 — j)/2(2 — j e ) t transforms Eqs. (fT2"j) 



and into Eq. (J22J) for 3A = 2A/a% and 7 = 7 e after fixing q = ^2(2 - %)/3(2 - 7). Its 
implicit solutions are obtained making these changes of variables in the solutions (|23 [) .(|24 [) 
and (f27)|) . (|2Tj) of the Sec. Ill A.. Explicit solutions also can be found for j e = 0, 1, 2 following 
the same steps as were done in that subsection. For a vanishing transverse pressure £ we 
reobtain Eq. (j22J) for the isotropic perfect fluid. 

To look into anisotropic fluid effects on the cosmological model, the behavior of the scale 
factors ()49|) will be investigated in two asymptotic regimes assuming that < 7 < 2 and 
7 e < 2. In the first regime, when a 3 < |cTq/2A| 1//( ^ 2 ~' /e ^ the contribution of the anisotropic 
fluid can be neglected and Eqs. (|5T|) . (jl2)l give t)~T and m . Hence, the scale factors 
(Jinj) behave as 

a % « a i0 T^ +qi+ W^ . (52) 

In the second regime, which starts when a 3 > |ag/2A| 1 ^ 2 ~ 7e ) the fluid dominates and Eqs. 
flT2"|) . (joTj) have the approximate solutions a 3 « uT 2 ^ e and m ^exp[(7 e /z/(7 e — 2))T ( - le ~ 2 ^' Ye ] J 
where v = (7 e v / 3A/2) 2,/7e is a constant. Inserting them in Eq. (j49|) . we get the large time 
behavior of the scale factors 



a i0 



v 



i+ 2r 



Therefore, this cosmological model never reaches an isotropic stage. For 7 e > 2 the existence 
of the solutions is linked to a negative energy density and the average scale factor has a finite 
time span. 



C. Anisotropic stiff fluid 

In this last subsection we will investigate the particular case of an anisotropic stiff fluid 
with equation of state p = p, where p is the average pressure, and find the general exact 
solution for a selected set of transverse pressures. The average index (J4H|) of this fluid is 
7 = 2 and the general solution of Eq. (JH4"j) is a 3 = a^t, where a is a constant. So the average 
expansion rate is H = l/3t. This suggests analyzing a transverse pressure which depends 
on the average scale factor as 
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where the transverse vector T and the parameter n are constants. The shear vector (jlOjl 
becomes 



a 



ej g+ (2-n)aga3(«-D' 71 ^ ' 



(55) 



(7 



— I r+ 



ln<r 



n = 2, 



(56) 



and Eq. (JBJ) gives the relation p = a® /3a 6 — a 2 / 2 between the energy density of the anisotropic 
stiff fluid and the shear scalar. 

Integrating the shear vectors (JHHj) and (JHHj) . we get the general solution 



flio a o t 3 +l3i exp 



(2 - n) 2 al n 



n^2, 



(57) 



a i0 ao t a o 



ai = a i0 a t ™o , n = 2, (58) 

For n < 2 the scale factors begin to evolve as it were an isotropic stiff fluid, a» ~ t^ +9i , 
(see Eq. (IHOjl ). But at late time the anisotropic stiff fluid dominates and, in this regime, 
the behavior of the scale factors is quite different from the isotropic case (see Eq. ()30p). 
For n > 2, the model has a high degree of anisotropy at early time, where the anisotropic 
stiff fluid dominates. After that, the scale factors evolve into those given by Eq. ([30)1 . This 
model could be used to describe a universe that is anisotropic at the very beginning. For 
n = 2, the solution behaves as solution (|3T)j) near the characteristic time t c = l/a . 



V. CONCLUSIONS 



We have found in the above results that the ISG of the Bianchi type-I cosmology is given 
internally by the Einstein gravitational action rather than being an additional external 
hypothesis imposed upon the model. Such a possibility appears to be reasonable for any 
unified gauge theory so that (a) the interactions are completely determined by the gauge 
invariance, and (b) a priori, the theory possesses an internal gauge group. 
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The relation between this ISG and the geometrical properties of the model is made 
through parameters such as the expansion scalar and the shear scalar, this means that the 
invariants of the group are expressed in terms of these quantities. The decomposition of 
the three-dimensional internal space as a direct sum of one-dimensional space containing 
the rotation axis and a two-dimensional orthogonal space allows one to write the Einstein 
equations in a simple manner. The three scale factors generate a vector representation of 
the ISG, and the associated transformations include a set that has been used to reduce 
the Bianchi type-I cosmology with a perfect fluid to a locally rotationally symmetric one 
with the same source and shear scalar. The existence of this particular representation is 
guaranteed by the existence of the internal space itself. 

We have given an integral representation of the general solution to the Bianchi type-I 
cosmological model with a perfect fluid and constant baryotropic index, and shown that the 
scale factors approach the average scale factor asymptotically for large times when < 7 < 2. 
These solutions begin as Kasner solutions, but after some characteristic time they have a 
final Friedmann-Robertson- Walker stage. Hence, the perfect fluid source dissipates the initial 
anisotropy of the model. 

In the last section we considered a Bianchi type-I cosmological model with three different 
types of anisotropic fluid; except for the anisotropic stiff fluid, the solutions found satisfy 
the constraint q ■ £ = 0. For an anisotropic fluid with constant stress-energy tensor we 
found a set of exact solutions which generalize those with a cosmological constant. When 
the magnitude of the transverse pressure is significantly smaller than the energy density of 
the fluid, the equation of state has a certain resemblance to that of an exotic fluid [20], 
called a modified Chaplygin gas equation of state. In the case where each pressure Pi of 
the anisotropic fluid is proportional to the energy density p, we showed that the transverse 
pressure increases the effective index and decelerates the average expansion of the universe. 
Finally, we found a general solution for an anisotropic stiff fluid with transverse pressure 
X = T/a 3n . For n > 2, the solutions are highly anisotropic at early time, but at late times 
they evolve into the isotropic stiff fluid solutions. 
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